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Abstract 

It has been shown recently that there exist analytical solutions of the Yang-Mills 
equations for non-Abelian Weizsacker- Williams fields which describe the distribution 
of gluons in large nuclei at small x. These solutions however depend on the color 
charge distribution at large rapidities. We here construct a model of the color charge 
distribution of partons in the fragmentation region and use it to compute the boost 
covariant momentum distributions of wee gluons. The phenomenological applications 
of our results are discussed. 



1 Introduction 



Recently, we constructed a QCD based model for the distribution of gluons in large 
nuclei at small x JE]. These fields are the non-Abelian analog of the well known 
Weizsacker-Williams fields in QED. A path integral was written down for the wee 
partons in the presence of valence quark sources. These sources are static charges on 
the light cone. In subsequent papers the saddle point solution of the path integral 
and the influence of quantum fluctuations about the saddle point solution were in- 
vestigated 0. It was shown later by one of us in Ref. (with J. Jalilian-Marian, 
A. Kovner and H. Weigert) that the validity of the original model could be extended 
to finite nuclei and even protons at very small values of x. Remarkably, it was also 
shown by J. Jalilian-Marian et al.|2|] and by Kovchegov || that the classical Yang- 
Mills equations for the non-Abelian Weizsacker-Williams fields can be solved and an 
analytical solution obtained for the nuclear gluon distribution at small values of x. 

However, the gluon distributions at a particular value of x depend on the color 
charge squared per unit area, x(y, Q 2 )i which is the integral of the color charge squared 
at all rapidities greater than the rapidity of interest. Formally, we define x t° be 



where fi 2 (r]',Q 2 ) is the color charge squared per unit area per unit rapidity at the 
rapidity r)' and Q 2 » Aq CD is the momentum transfer squared of an external probe. 
Note also that the space-time rapidity is defined by the relation 



Here x R ~ 2R/ \fl r ) is the Lorentz contracted thickness of the nucleus in the infinite 
momentum frame (7 >> 1) and x~ = (t — z)/\/2 is the longitudinal light cone co- 




(1) 



r) = r] R + log(x R /x ) . 



(2) 
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ordinate. At central rapidities, the space-time rapidity is approximately equal to the 
momentum space rapidity. 

In order to explicitly compute the parton distributions in nuclei at small x, we 
need to model the color charge distribution of partons in the fragmentation region. 
In this paper we will construct a simple model which captures the essential physics 
of the fragmentation region in nuclei. The expression for \ that results will lead 
to a boost covariant expression for the gluon distribution in large nuclei at small 
x. A conceptual issue that is resolved as a consequence is the description of the 
non-Abelian Weizsacker-Williams fields in the nuclear rest frame. This enables us in 
principle to compare our results to calculations in the nuclear rest frame ||. 

In large part, this work is an extension of the nice analysis done by Kovchegov || 
who, in a simplified model of a nucleus, formulated in the nuclear rest frame, argued 
that the nucleus acts as a source of color charge. Here we incorporate a model of 
the longitudinal dynamics and construct a solution in the fragmentation region of the 
nucleus. We argue that in most of the fragmentation region for a sufficiently large 
nucleus, our weak coupling methods are valid. 

We should warn the reader, that our results for realistically sized nuclei are semi- 
quantitative and qualitative at best. For realistically sized nuclei, the typical energy 
scale in the problem is of the order of 300-400 MeV at RHIC energies and 750-900 
MeV at LHC energies at central rapidities ||, which is marginally large enough (for 
our weak coupling calculations to be valid) only in the case of the latter. In order to 
get a scale well in the weak coupling regime would require a nuclear baryon number a 
couple of orders of magnitude larger than that of the largest nucleus or equivalently, 
much higher energies. Nevertheless, we believe the insight we gather from this picture 
may be useful for obtaining a reasonable, QCD motivated conceptual framework in 
which to think about high energy nuclear collisions. 
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Our results also enable us to study the sytematics of nuclear gluon shadowing 
and other phenomenological observables at low x. These will be especially relevant 
for the deep inelastic scattering experiments off nuclei which have been proposed at 
HERA ||. A caveat is that our result for the gluon distribution is classical. However, 
another important result of Jalilian-Marian et al. || is that the effective action for 
wee partons demonstrates a self-similar structure as we go to lower rapidities. The 
effective color charge x obeys a renormalization group equation analogous to the 
DGLAP or BFKL equations [§] (we note that in a recent paper, it has been shown 
explicitly that the BFKL kernel is recovered in this approach [[Hj). Therefore, at 
each rapidity slice, the gluon distribution is given by the classical gluon distribution 
with the quantum effects absorbed in the effective scale \ at that rapidity slice. The 
residual quantum corrections are of order as which is much smaller than unity for 
the rapidity slice of width Ar]. These can in principle be computed. In other words, 
the gluon distributions computed using our model for x, provide the initial conditions 
for the renormalization group equations which can be solved self-consistently for \- 

The rest of the paper is organized as follows. In section 2, we briefly discuss and 
summarize the results of Ref. || for the classical gluon distribution. In section 3, 
we construct our model for x, the color charge squared per unit area. Our model 
depends on one free parameter, k, which is the average color charge squared per 
nucleon integrated over one unit of rapidity. Section 4 contains a discussion of how 
k may be determined and the implications of our model for shadowing in very large 
nuclei. The final section contains a brief summary and points out directions for future 
work. 
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2 An analytical solution for correlation functions 
of non— Abelian Weizsacker— Williams fields 



In Ref. [[J £5|, a path integral was written down for wee parton modes in a nu- 
cleus with x << A~ x l z and transverse momentum Aq C d « h « y/x- Here 



X = dr]'fi 2 (rj', Q 2 ), where /i 2 is the average color charge squared per unit area at a 
particular rapidity rf and a transverse momentum resolution Q 2 . The latter is nor- 
mally associated with the momentum resolution of an external probe. The measure 
for these modes (in light cone gauge A + = 0) is 



Above, y has the same definition as in Eq. 0. In the trace that defines the path integral, 

the sum over the valence quark color charges has been replaced by a Gaussian measure 

with the weight p 2 . The action Sqcd above refers only to the pure gauge degrees of 

freedom. The above action does not include sea quarks (which are as suppressed) 

but is easily extended to do so. Since /i 2 is the only scale in the problem, the coupling 

constant will run as a function of this scale. For p 2 » A.q CD , our effective theory 

will be in the regime of weak coupling. 

This procedure is justified when the following two assumptions apply Q The first 

is that the valence quarks are static charges on the light cone with J +,a = p a (xt)5(x~) 

as the only large component of the valence quark current. All other components of 

the current are suppressed by 1/ ' P + . The second assumption is that the wee partons 
1 For a recent discussion of a gauge invariant form for the above action, see Ref. [[Io|. This 
discussion while of general interest, is not applicable here since we restrict ourselves to the semi- 
classical result. 



z 




(3) 
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couple to a large color charge which is classical. What this implies is that the charge 
squared per unit area, < Q 2 >, that the wee partons couple to is given by the relation 
< Q 2 >= N g < Q 2 > qU ark, where N q is the number of quarks per unit area (oc A 1 / 3 ) 
and < Q 2 > qua rk is the charge squared per unit area due to a single quark. It can 
be shown that the central limit theorem can be applied to this problem and that the 
corrections to the classical result (which may in principle be computed) are of the 
order 1/ \jN~ q ~ A -1 / 6 for the root mean square classical charge |J. For an alternative 
derivation of the Gaussian measure in a specific model, see Ref. || . 
To compute the classical nuclear gluon distribution function 
(j at o r 

~¥k = J2^ 2lk+l J Sxd " x ' eiHx ~ X,) < A a t {x-,x t )A){x'-,x' t ) > p , (4) 

one needs to compute the saddle point solution of Eq. [|, or equivalently, solve the 
classical Yang-Mills equations for the source J + defined above and finally average 
over the product of the classical fields at two space-time points with the Gaussian 
measure in Eq. ||]. Note that the brackets < • • • > p above correspond to the Gaussian 
averaging over the static, light cone sources. One finds a particular solution where 
A ± = and AHi = 1,2) is a pure gauge field which satisfies the equation 

Di^j- = gp(v,x±) , (5) 

where Di is the covariant derivative [I]. 

In Ref. 0, the above equations were solved to obtain A % as a function of the color 
charge density p and the distribution function computed. Here we will merely state 
the result of Ref. [0] and refer the reader to that paper for the technical details. The 
distribution function for the non-Abelian Weizsacker-Williams fields in coordinate 
space is 

dN 4 (iV c 2 -l) / , 2a2 ^ x( ,,q> 



d 2 x t drj g 2 N c xf 



1 - {^lco)™ mV " (6) 



The momentum space distributions are easily obtained by Fourier transforming 
Eq. §. The result is 

1 dN 4 (N 2 — 1) 1 /-iog(V^) 

x / -0± J (k t u ± ) 1 - (u ± A QCD f] . (7) 

Jo Uj_ 1 J 

Above, r = 1.12 fm, is the nucleon mass and (3 = 2ira 2 ^N c x{y',Q 2 )u 2 ± . For 

as^/x << the distribution has the familiar Weizsacker- Williams form oc g 2 /k 2 , 

while for Aq C d « h « ois^fx, the distribution is oc 4y \og(k 2 / g 4 N c x(y)) ■ Even for 

rapidities far from the beam rapidity, it is reasonable to assume that y spa cetime ~ Vmom 

for wee glue and light sea quarks. For kx very small compared to Ois^fx, the above 

formula may break down, but we expect that a large number of physical quantities 

will not be infrared sensitive. 

While the above expressions provide a very interesting picture of the behaviour of 

the gluon distribution as a function of the transverse distance Xt, the dependence of 

the distributions on the spacetime rapidity is hidden in the function x(Vi Q 2 ) which 

has not been specified. In order to specify x f° r ai l rapidities r], we need a model of the 

nuclear color charge distribution (at large Q 2 ) in the fragmentation region as well. In 

the next section, we construct a model for x at some reference scale Ql >> Aqcd- As 

we will discuss below, x obeys evolution equations in both y and Q 2 . The important 

point to note here is that with our model of x the expression in Eq. |6] is then boost 

covariant. This enables us to describe what the wee parton distributions look like in 

any frame. 

Finally, the parton distributions computed in the above described manner can 
be related to the empirically extracted gluon distributions at the scale Q 2 by the 
equation 



This simple relation is obtained using light cone quantization and light cone gauge 
A+ = 0|. 



3 The rapidity distribution of color charge in nu- 
clei at a fixed resolution scale Qq 

In this section we shall formulate a model for x, the color charge squared per unit 
area per unit rapidity, integrated upto the rapidity of interest. A similar model for x, 
albeit without the explicit y dependence, was discussed recently by one of us and M. 
Gyulassy |J. A detailed model in the fragmentation region is relevant for studying 
the renormalization group evolution of x towards smaller rapidities and larger Q 2 . 

Consider a nucleus at rest. For a large nucleus, the color charge squared per unit 
longitudinal distance z, dC/dzd 2 x± is a constant within the nucleus when measured 
on longitudinal resolution scales which are large compared to the inter-nucleon sep- 
aration. The color charge squared is the typical value of charge fluctuations in the 
fragmentation region of the proton-about 1 unit of rapidity. Here we must include 
the fluctuations due to the gluons as well as the quarks, so this may be roughly a 
factor of two larger than is the case for just the valence quarks. Note also that we get 
a non zero color charge only because we are choosing to measure the charge density 
on a transverse size scale much less than a fermi. In a large nucleus, the color charge 
resolved by the probe corresponds to a higher dimensional representation of the color 
algebra and is therefore classical JI|, [|. The probe sees only small transverse bits of 
each nucleon but the fluctuation fields add incoherently from nucleon to nucleon. 

The non-Abelian Weizsacker-Williams fields appear as solutions in the nuclear 
rest frame to an external probe in the following sense. Outside each nucleon, for the 
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long wavelength small x modes, there is only the gauge transform of the vacuum field. 
As we pass from one side of a nucleon to the other, we get a different gauge transform 
of the vacuum. The field inside a nucleon is not known. Further, since the field is 
pure gauge outside each nucleon, or in other words, the field strengths are zero, the 
classical force exerted by the long wavelength modes of one nucleon on a neighboring 
one is negligible. We therefore have a model of the nucleus with short range nuclear 
forces where it is the short wavelength (x > 0.1, where x ~ 0.1 is the wavelength 
corresponding to the inter-nucleon separation) modes which are responsible for intra- 
nuclear dynamics such as nuclear binding. This model is therefore consistent with 
the conventional picture that nucleonic effects in nuclei (such as the EMC effect 
for instance) are unaffected by the physics of small x modes in the nucleon (for a 
discussion see for instance Ref. ||12||). The typical change in the Weizsacker- Williams 
field seen by the external probe should then be described by a stochastic source with 
a strength which must be determined phenomenologically by an empirical knowledge 
of the nucleon structure functions. It is with this idea in mind that we formulate the 
model below. 

In the nuclear rest frame we choose to localize the color charge as 

An 

^oc6(z + 2R)6(-z), (9) 
dz 

where C is the classical color charge squared per unit transverse area at the scale Qq 
of the external probe. In the boosted frame, recall that we define r] = r]R + ln{x^ l /x~) 
where x~ = (t — z)j\[2 and x^ ps 2R/\/2^. 

It is easily seen that the above charge distribution is independent of the light cone 
time x + for x + << 7 and 7 >> 1. This result is obtained by boosting the field in 
the rest frame to the infinite momentum frame. In the classical limit we are working 
in, the field in one frame is just the Lorentz transform of the field in another. In this 
limit, the commutators of the color charges vanish, thereby allowing us to relate the 
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charge distributions in the two frames. 

The color charge squared per unit transverse area per unit rapidity in the nuclear 
rest frame is 

dC dC dz 
drj dz drj ' 

and from our definition in Eq. we have 
for < rj < oo. 

The color charge squared density in the nuclear rest frame is 

dC . , 

- = Kpo (12) 

where k is the typical charge squared per nucleon in one unit of rapidity and po is 
the density of nucleons. The latter may be chosen for a large nucleus to have the 
typical value for nuclear matter density. All of the uncertainty about the color charge 
distribution inside each nucleon is therefore represented in k. 

We can now boost all this to the frame where the nucleus moves very fast. The 
main effect of the boost is to replace the density po by the nuclear gamma factor times 
Po- We find therefore that the new charge squared per unit area per unit rapidity 
(which is the p 2 introduced in the previous section) becomes 

II 2 (V,QD = 2RK Po e^ ;v>VR 

= ',7}<7] R . (13) 

This rapidity distribution is amusing. In the rapidity variable, most of the color 
charge squared is concentrated in one unit of rapidity above the beam rapidity. This 
is a boost invariant statement. The region where the charge squared distribution falls 
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to that typical of a single nucleoli is however Arj ~ ln(R/l Fm). Even though most 
of the nucleons are concentrated in one unit of rapidity, it takes the order of ln(R) 
units in rapidity before the charge squared has fallen to a small value! 

It is also interesting to consider how this picture of the fragmentation region can be 
formulated in terms of the momentum space rapidity. A similar argument to the one 
above applies if we formulate our discussion in terms of a momentum space rapidity 
that we define as y = yn + \og(k + /-Pr). Note that this definition is not the usual one 
since P^ is not the nuclear momentum per nucleon P + , but instead the momentum 
that is conjugate to the Lorentz contracted width x^ of the nucleus. This definition 
has the following advantages. Firstly, from the uncertainty principle, it is equivalent 
to the space-time rapidity we defined previously since log(k + /P^) = \og{x^/x^). 
Secondly, by this definition, momenta P^ < k + < P + correspond to fragmentation 
region rapidities, or in terms of x~, wavelengths shorter than the Lorentz contracted 
width of the nucleus. One therefore expects that the charge distribution in terms 
of the above defined momentum space rapidity is identical to the distribution in 
Eq. [L3|: xijliQ 2 ) = xhjiQ 2 )- What fraction xr of the nuclear momentum P + 
does P£ correspond to? It may be determined by requiring that it correspond to a 
wavelength equal to the Lorentz contracted width of the nucleus = x^. Hence, 



where is the nucleon mass and r = 1.12 fm. For A = 200, this corresponds to 
xr ~ 0.025. Therefore, by our definition, for A = 200, fragmentation region rapidites 
y > yji correspond to light cone momentum fractions x > 0.025. Of course even for 
these large nuclei our classical considerations are very marginal and can be safely 
applied in the theoretical limit of much larger nuclei or equivalently, when the parton 
densities in the fragmentation region are already large. 

Returning to Eq. |13|, because most of the color charge squared is contained in one 




(14) 
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unit of rapidity, there is not much dynamical evolution of the glue distribution due to 
the glue from different nuclei interacting in the fragmentation region. This is because 
the typical interaction strength is on the order of a$ which is small due to the high 
field density and because the typical distance it takes to evolve the gluon distribution 
due to the interactions of the gluons among themselves is of order Treating 
the nucleons as a source of glue in the nuclear fragmentation region is therefore 
justified so long as sources due to gluons of order one unit of rapidity from the beam 
rapidity are added to those from the valence quarks. This of course introduces a 
phenomenological parameter (k above) which we cannot precisely determine from 
first principles. However, as will be discussed below k can be determined empirically 
from the nucleon gluon distributions. 

The gluon rapidity distribution varies slowly (oc as per unit rapidity) when 
Ay << l/a,s away from the beam rapidity. It is obtained by solving the classi- 
cal equations of motion for the source distribution above. Integrating the expression 
in Eq. [L3| above, we obtain 

X(y, Ql) = 2 Rp K ;y<y R 

= 2Rp Ke yR - y ;y>y R , (15) 

where, we remind the reader, y R is the momentum space rapidity defined previously. 
The ratio of the nuclear gluon distribution to the nucleon parton distribution times 
A in the nuclear fragmentation region for very large nuclei is plotted in Figure 1. Our 
classical considerations result therefore in a charge distribution that is independent 
of rapidity for y < y R . Since any gauge invariant configuration of fields must satisfy 
Gauss's law, the contribution of glue at small x is necessarily associated with the 
valence degrees of freedom and therefore also exhibits boost invariance. This is also 
clear from examining Eq. ^. 
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4 Discussion 



From our elementary considerations above, we have obtained an analytic, boost co- 
variant form for the color charge squared distribution \ in the nuclear fragmentation 
region at the momentum resolutution scale Q\. It relies on one parameter k which 
encapsulates our uncertainty about the non-perturbative physics in the fragmenta- 
tion region. However, since k has the transparent physical interpretation of being 
the color charge squared per nucleon, some of the uncertainty can be eliminated by 
expressing k in terms of the nucleon gluon density. Hence, at some reference scale 

k = ^- dx (V(x, Ql) + S(x, Ql)) + j^-j dxG(x, Ql) . (16) 

Here V, S and G denote respectively the valence, sea quark and gluon densities. They 
are integrated over the roughly one unit of rapidity that one estimates the nucleons to 
occupy in the fragmentation region. The factors multiplying the parentheses above are 
the color charge squared per quark and the color charge squared per gluon respectively. 
The last term for example, has the interpretation that it is the color charge squared 
per gluon times the number of gluons in the nucleon contained in one unit of rapidity 
in the fragmentation region. Similarly so for the first term. 

The parton densities at the reference scale Qq can be parametrized by fits to the 
nucleon deep inelastic scattering data. A convenient one is the NLO parametrization 



of Gliick, Reya and Vogt Jy|. In principle, the sensitivity of the results to different 
parametrizations can be investigated. As is well known by now, the GRV parametriza- 
tions of the data at the (extremely low!) reference scale {Q g $ v ) 2 = 0.34 GeV 2 , evolved 
with the next to leading order Altarelli-Parisi evolution equations, agree well with 
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the DIS data from HERA. The GRV parametrizations have the generic form 

xf(x, Q 2 ) = Nx a (l + Ax b + Bx + Cx 3/2 ) (1 - x) D , (17) 

where f = V, S,G and the coefficients and powers in the above expression are polyno- 
mials of log(log(Q 2 / Aq CD )) . The detailed expressions are given in Ref. flT5| |. For our 



classical considerations to apply, we require the reference scale Q » Qq TV ~ h-QCD- 
The parton densities which determine k can then be parametrized by the GRV parton 
distribution evolved to the scale of interest Qq. 

With the initial data for x as input, the Wilson renormahzation arguments of 
Ref. can be used to determine it at higher Q 2 and smaller values of y. It has the 
form 

d X N c 1 dN 



dydQ 2 N 2 - 1 ttR 2 dydQ 2 ' v ' 

where dN/dydQ 2 is the distribution function in Eq. [7]. This renormahzation group 
equation can be formulated in DGLAP form by first integrating over y, 

dx N c 1 f °° dN 

d V 17717^ ■ ( 19 ) 



dQ 2 N 2 - 1 ttR 2 J y * dy'dQ 2 
Integrating over Q 2 , we obtain the BFKL-like form 

d -* = J^^{ Q2 dQ> 2 ^L- (20) 
dy N 2 -l7rR 2 J ^ dydQ' 2 ' 1 } 

It was shown recently that for large momenta the kernel of the equation satisfied 
by the right hand side of the above is precisely the BFKL kernel [PH . However in 
general, as discussed earlier, the gluon distribution on the right hand side of the above 
equations has a non-linear dependence on x- 

A clear picture of shadowing then emerges in the above model. From the pertur- 
bative QCD standpoint, one of the possible sources of shadowing is from the the input 
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parton distributions which are modified upon Q 2 evolution by the splitting functions 
a la Altarelli-Parisi [16], ^Ofl . The other source is due to recombination and screen- 



ing effects which enter into the evolution equations as higher twist effects [17, 18 



There have been a number treatments of shadowing with take both these effects into 



account [19|, g0| |21| . It is often believed that current data can be accounted for by 
shadowing of the initial (non-perturbative) parton distributions and leading twist 
DGLAP evolution of these distributions. 

On the other hand, our model in the previous section suggests that for sufficiently 
large nuclei i) the initial parton distributions can be determined perturbatively (mod- 
ulo k) and ii) they are not shadowed. However the shadowing away from the frag- 
mentation region is large due to higher twist effects which, in this model, are already 
present at the classical level. If we substitute the classical result of Eq. [F] into the 



right hand side of Eq. 19 



I* = ik r £ t mq z) 11 - exp < 3 • <2i) 

where (3 = 2ira 2 ^N c x{x' , Q 2 )z 2 . If we expand out the 1 — exp(- ■ •) term on the right 
hand side and make the identification x — jfftfxG(x, Q 2 )/tcR 2 , we obtain an evolu- 
tion equation in the gluon density G(x, Q 2 ) which on the right hand side includes all 
powers of G(x,Q 2 )/Q 2 R 2 -oi all twists. The first term in the expansion is just the 
usual DGLAP evolution equation at small x. A very similar expression to the above 
has been obtained by Levin and collaborators working in the nuclear rest frame-and 
is the perturbative QCD form of the Glauber-Gribov-Mueller formula of high en- 
ergy scattering ||22|| . Since the second term in the expansion above can be related 
to diffractive scattering, diffractive scattering data can be used to determine the de- 
gree of nuclear shadowing in this approach. This was done recently by Capella et 
al. who used HERA diffractive data on e-p scattering to reproduce shadowing data 
from fixed target experiments-thereby providing a nice check of the Glauber-Gribov- 
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Mueller expression upto second order in the gluon density p3| . 

Ultimately, when parton densities are large, the quantum corrections to the above 
formula will be large and have to be computed. For instance, the second term in the 
above equation gives only the recombination (or diffractive) contribution to shadow- 
ing (with positive sign!). The screening corrections will appear when one includes 
quantum corrections. If the well known AGK sum rules |24j] apply, then as argued 
by Mueller and Qiu [18|]- the expression will be the same as the above but with a 
negative sign! This has to be confirmed in our approach before anything quantita- 
tively new can be added to current phenomenology. Moreover, we hope that in this 



effective action approach it will be easier to compute higher order terms as well [C 



5 Summary and Outlook 

In this paper we have formulated a simple model for x{Hi Q"o)i ^ ne c °l° r charge squared 
integrated from the rapidity of interest upto the beam rapidity at some reference 
scale Qq. At this scale, x is the color charge in the nuclear fragmentation region. 
X depends on the non- perturbative physics at this scale through k which is the 
charge squared per nucleon and is simply related to the nucleon parton distributions. 
The expression derived for x is boost covariant. The Wilson renormalization group 
equation discussed in Ref. can be used to evolve x to Q 2 > Ql and/or rapidities 
away from the fragmentation region. We have argued that this classical picture 
provides important insight into the nature of perturbative gluon shadowing and have 
argued in the previous section that our results are closely related to similar approaches 
in the nuclear rest frame. This is indeed as it should be because the model we 
have constructed is boost covariant. The boost covariant gluon distributions thereby 
constructed can be applied to a wide range of phenomenology including gluon and 
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quark shadowing [II, [13[, intrinsic quark production |^6| and nuclear collisions at 
ultrarelativistic energies |27, |5|. It is however important to understand the quantum 
corrections in the infinite momentum frame effective action picture before significant 
additions can be made to current phenomenology p5|. 
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6 Figure caption 

Fig. 1: Ratio R f ™ 9 = xG a {x,QI)/Ak (where k is defined as in Eq. |16|) CIS db function 
of rapidity in the fragmentation region. Above is defined by Eq. |14|; ybeam — Ur ~ 

log^/S) and Q2 A 2 Qc ^ 
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